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ABSTRACT 

The prototype of accreting, pulsating white dwarfs (GW Lib) underwent a 
large amplitude dwarf nova outburst in 2007. We used ultraviolet data from 
GALEX and ground-based optical photometry and spectroscopy to follow GW 
Lib for three years following this outburst. Several variations are apparent during 
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this interval. The optical shows a superhump modulation in the months following 
outburst while a 19 min quasi-periodic modulation lasting for several months is 
apparent in the year after outburst. A long timescale (about 4 hr) modulation 
first appears in the UV a year after outburst and increases in amplitude in the 
following years. This variation also appears in the optical 2 years after outburst 
but is not in phase with the UV. The pre-outburst pulsations are not yet visible 
after 3 years, likely indicating the white dwarf has not returned to its quiescent 
state. 

Subject headings: cataclysmic variables — stars: individual (GW Lib) — stars: 
oscillations 



1. Introduction 

GW Librae was first thought to be a nova (Duerbeck 1987) because of the large amph- 
tudc of an outburst during its discovery in 1983, but a later spectrum at quiescence revealed 
it to be a dwarf nova. The lack of outbursts for the next two decades, combined with the 
large outburst amplitude argued for a classification as a WZ Sge type of dwarf nova (How- 
ell, Szkody & Cannizzo 1995). Time-resolved spectroscopy (Szkody, Desai & Hoard 2000; 
Thorstensen et al. 2002) revealed a very short orbital period of 76.78 min, consistent with 
this classification. However, the unique nature of GW Lib became apparent when it was 
discovered to have non-radial pulsations (Warner & van Zyl 1998) and it became the pro- 
totype of the small class of accreting, pulsating white dwarfs which currently contains 13 
members. Exploration of its amplitude spectrum by van Zyl et al. (2000, 2004) identified 
three prominent pulsation periods at 648, 376 and 236 s along with fine structure of closely 
spaced frequencies around the 648 and 376 s periods. The amplitudes of the pulsations 
changed on monthly timescales, typical of most hydrogen atmosphere white dwarf pulsators 
(ZZ Ceti stars). Ultraviolet observations of GW Lib (Szkody et al. 2002) revealed a hot 
white dwarf with a temperature consistent with the theoretical instability strip at 15,000K 
due to He II ionization (Arras et al. 2006). The best fit was achieved with a high mass 
(0.8 Mq) white dwarf with 0.1 solar metal abundance and a dual temperature (63% of the 



^ Based on observations made with the NASA Galaxy Evolution Explorer and with the Apache Point 
Observatory (APO) 3.5m telescope. GALEX is operated for NASA by the California Institute of Technology 
under NASA contract NAS5-98034. APO is owned and operated by the Astrophysical Research Consortium 
(ARC). 
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surface at 13,300K and 37% at 17,100K). The UV data showed higher amphtude pulsations 
than at optical wavelengths with the same periods as evident in the optical. All stellar pul- 
sators, including variable white dwarfs, reveal wavelength dependent amplitudes whenever 
their atmospheres suffer from limb-darkening effects. 

Recent studies of the temperatures of nine of the other accreting pulsators (Szkody et 
al. 2010) showed that their instability strip is wider than for ZZ Ceti stars and some stars 
appeared to stop pulsating during the Hubble Space Telescope observations. One possible 
explanation for at least one of the systems that did not show pulsations was that accretion 
heating from an outburst could have caused the white dwarf to move out of its instability 
strip. Outbursts are known to cause major heating of the white dwarfs (e.g. WZ Sge; Godon 
et al. 2006) which then take several years to cool to their quiescent values, even though the 
optical light reaches quiescence in only a few months. 

On 12 April 2007, GW Lib went into outburst (Templeton et al. 2007), reaching 8th 
magnitude from its quiescent brightness near 17th magnitude. The outburst was followed 
spectroscopically (Nogami et al. 2009; Hiroi et al. 2009; van Spaandonk et al. 2010a) and 
photometrically (Copperwheat et al. 2009, Schwieterman et al. 2010) as well as with X-ray 
telescopes (Byckling et al. 2009). The spectroscopy revealed the usual prominent optically 
thick disk near outburst with a slow return toward quiescence over 3 months. However, the 
outburst was unusual compared to other dwarf novae and WZ Sge objects in several ways. 
The X-ray flux of GW Lib increased by 3 orders of magnitude during the outburst and 
remained one order of magnitude higher than pre-outburst even at 2 years past outburst. As 
Bykling et al. (2009) point out, GW Lib is the only other dwarf nova other than U Gem that 
shows a larger X-ray flux during outburst than quiescence. The optical flux of GW Lib also 
remained about a half-magnitude higher than quiescence 2-3 years past outburst, whereas 
most WZ Sge stars return to optical quiescence in a few months (e.g. AL Com and WZ Sge; 
Szkody et al. 1998). Finally, there was no evidence of the echo outbursts that many WZ 
Sge stars show on the dechne from outburst (e.g. EG Cnc; Patterson et al. 2002). 

We observed GW Lib for three years following its outburst, using ultraviolet data from 
GALEX, and ground-based optical photometry and spectroscopy. Our aim was to under- 
stand the heating and cooling of GW Lib as it moved out of its instability strip and then 
re-entered. As of the current time, the white dwarf has not yet resumed its pre-outburst 
character. Yet, the photometry has revealed some interesting facets of GW Lib during its 
decline from superout burst. 
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2. Observations 
2.1. GALEX Data 

Ultraviolet data were obtained using the GALEX satellite in April- June of 2007, 2008, 
2009 and 2010 (Table 1 lists the times). Both the near-ultraviolet (NUV: 1750-2800A) and 
the far-ultraviolet (FUV; 1350-1750A) detectors (Martin et al. 2005) were used, except for 
the 2010 data when the FUV detector was no longer operational. The mean magnitudes from 
each set of observations were used to construct a light curve over the three years (Figure 1). 
From the GALEX online documentatioij^, the zeropoint of the AB magnitudes is given by 
mo= 18.82=1.40x10"^^ ergs cm'^ s"^ A"^ for FUV and mo= 20.08=2.06 xlO'^^ ergs cm'^ 
s-i A-i for NUV. 

In order to be able to analyze the data for any short-term periodicity related to pulsa- 
tions, the time-tag data were binned (into 29s exposures for the 2007 data and 59s for the 
2008-2010 data), calibrated and measured with the IRAI^ qphot routine. An 8 pixel radius 
aperture was used for the stars and an annulus of 10-13 pixels for sky measurements. Nine 
reference stars in the field were measured in the same way as GW Lib to serve as a check on 
any variations that could be due to the various satellite orbits. 

2.2. Optical Photometry 

Optical photometric data were obtained with multiple telescopes between 2007-2010, 
beginning four days after the ouburst of GW Lib in April (Templeton et al. 2007). The 
primary telescope used was the robotic 1.0 meter of New Mexico State University (NMSU) 
which utilizes an E2V 2048x2048 CCD. In 2007, the night typically began with 2-6 observa- 
tions through U,B,V, R and / filters. Longer observing runs were then obtained through the 
U filter. Beginning in 2008, a BG40 filtecl was also used and in 2009, only the BG40 filter 
was employed. This filter is a Schott glass filter with a broadband response from 3500-6000A 
which allows better S/N than the U filter, while still permitting better detection of pulsa- 
tions which have higher amplitude in blue light than in red (Robinson et al. 1995). Table 



^http:galexgi. gsfc.nasa.gov/tools/index. html 

■^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

''http://www.us.schott.com/advanced_optics/enghsh/download/schott_bandpass_bg40_2008_e.pdf 
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2 lists the ranges of exposures for each night, with the U band images occupying the upper 
end of the range. Throughout each night, the exposure times were automatically adjusted 
to compensate for changing conditions such as airmass and cloud cover. Differential pho- 
tometry with respect to calibrated standard stars in the field was used to obtain magnitudes 
with resulting photometric errors of up to a few hundredths of a magntitude. 

On 2007 May 29, a datasct was taken with the 3.5m Wisconsin, Indiana, Yale, National 
Optical Astronomy Observatory (WIYN) telescope through a BG39 filter which has similar 
transmission to the BG40. Observations were made with the MiniMo CCD with 2x2 binning 
and differential photometry was used for the reductions. 

Observations with the 3.5m telescope at Apache Point Observatory (APO) were made 
on 2007 June 1 and 2008 March 28. The data were taken through a BG40 filter with the 
Agile CCD and differential photometry was again used to obtain the magnitudes. 

Data from the Im telescope at the United States Naval Observatory Flagstaff Station 
(NOFS) were also obtained in 2009. A high speed camera was used on May 15 with a 
V filter, and the Tektronix 1024x1024 CCD with B filter was used on May 16 and 17. 
Additionally, the Sonoita Research Observatory (SRO) 0.5m telescope was used on 2009 
May 16. With this telescope, unfiltered images were taken with an SBIG STL 6303 CCD. 
Differential photometry was employed in both reductions. 

Southern hemisphere observations were made with the 4.1m Southern Astrophysics Re- 
search (SOAR) telescope on 2009 May 17. A combination of two E2V CCDs were used to 
make measurments through a V filter. Differential photometry using the IRAF qphot routine 
was used to obtain magnitudes. A second southern dataset was obtained with the South 
African Astronomical Observatory (SAAO) 1.0 m Elizabeth telescope beginning 2009 May 
18. A B filter was used with the STE4 CCD camera. Differential photometry was then used 
to obtain a light curve. 

In 2010, two nights of photometry were obtained on March 4 and 5, using the KPNO 
2.1m telescope with a BG39 filter and the STA2 CCD. As with the other reductions, differ- 
ential photometry was used. 

Approximate calibrated B magnitudes were obtained for all nights through a variety of 
methods. Generally, B and V magnitudes for comparison stars in the field were obtained 
from the AAVSO Variable Star Database and offsets computed to then cahbrate the GW 
Lib lightcurves. With the NMSU data, offsets for each night could usually be obtained from 
the B filter. For the 2007 nights with no B filter, the offset was determined by comparison 
of the mean magnitude to the corresponding magnitude from the AAVSO lightcurve. On 
2008 Jun 20, data were taken in both B and BG40 filters, which enabled a determination of 
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the offset to transform the BG40 magnitudes of 2008-2010 to approximate B magnitudes. 
The AAVSO standards were also used to obtain a transformation of the clear filter of SRO 
to an approximate B magnitude. To tranform the V filter data of SOAR and NOFS, it was 
assumed that the B — V color of GW Lib was zero. 



2.3. Optical Spectroscopy 

An optical spectrum of GW Lib was obtained on 2008 April 17, using the R-C Spec- 
trograph and grating KPGLl-1 on the CTIO 4m telescope. Two 30 min exposures were 
medianed together. The slit width was 1 arcsec, comparable to the seeing (determined from 
a slice across the spectrum). The airmass at the time of observation was 1.04 and the blue 
DA white dwarf EG274 (Hamuy et al. 1992) was observed at airmass 1.03 immediately 
following GW Lib. Due to the low airmass and the closeness in time of the standard, we 
estimate the flux uncertainty to be about 10%. The resulting cahbrated flux at 5500A trans- 
lates into a visual magnitude of 16.3, which is within the range of the AAVSO light curves 
during 2008. 

Table 2 gives a summary of the dates and times for all the optical data. 

3. Light Curves 

Figure 1 shows the mean of the ultraviolet and optical blue light curves from the obser- 
vations. The optical light undergoes a rapid linear decline in the first month after outburst, 
followed by a steep drop of 3 magnitudes between May 7-10 (days 227-230 on the plot) 
and then a slower decline continuing for years after outburst. The steep drop is generally 
ascribed to large changes in the disk that bring it close to its pre-outburst state (Byckhng 
et al. 2009). The GALEX cover&ge began after this steep drop had occurred. The ultravi- 
olet brightness dechnes by 1.5 magnitudes from 2007 June to 2008 June, while the optical 
changes by less than a magnitude during this time (the low point at day 560 is likely due 
to using a redder standard for calibration due to the smaller field of view of the detector 
during that observation as well as a different filter). The difference between the GALEX 
FUV and NUV magnitudes decreases from 0.4 mag in 2007 to 0.08 mag in 2009. The optical 
blue magnitude decreases slightly in 2010 while the ultraviolet NUV magnitude remains the 
same. These dechnes are consistent with the coohng of a hot source that peaks in the FUV, 
while the NUV and optical comprise the tail of the flux distribution. The similarity of the 
NUV fluxes from 2008-2010 indicates the continued presence of a hot source. 
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Figure 2 shows the individual GALEX hght curves in NUV and FUV filters for the 
various satellite orbits on a given date, while Figure 3 shows the optical light curves for the 
nights with data coverage longer than an hour. 

The individual GALEX hght curves are flat throughout the 2008 May and June obser- 
vations. The optical hght curves are also relatively flat from 2007 April 16-22 (Figure 3a), 
which corresponds to the timescale of the linear decline in Figure 1. The flat light curves are 
consistent with the domination of the system light by the optically thick uniform accretion 
disk. On May 7, at the start of the steep optical decline (Figure 1), the optical light curves go 
through major changes (Figure 3b). On May 7, the decline starts to appear during the 4 hrs 
of observation. Three nights later, at the end of the rapid decline on May 10, the light curve 
shows large fluctuations of several tenths of a magnitude, with a double-hump appearance 
repeating at a timescale near 80 min. On May 14, a slow sinusoidal type modulation is ap- 
parent near 80 min that is also visible on May 29. This is the same interval in which Kato et 
al. (2008) find a stable period of 77.98 min which they call a "late superhump" period. The 
erractic fluctuations in our May 10 data may be the disk transitioning from its early to late 
superhump state. Kato et al. (2008) note the strange stability and longer than usual period 
increase of this superhump compared to normal WZ Sge stars. Copperwheat et al. (2009) 
note a similar long period in their 2007 July data but with a lower amplitude modulation 
than we see in May. These longer period superhumps have been postulated to be due to 
some amount of gas being outside the 3:1 resonance radius, which causes a tidal instability 
leading to an eccentric disk (Whitehurst 1988). The lack of modulation of the NUV and 
FUV data during this time imphes that the inner, hotter disk areas are not affected. As the 
superhump amplitude decreases with time, the amount of mass in this outer area of the disk 
may be shrinking or there may be less irradiation of the disk as time progresses, although 
the low inclination of GW Lib (~ lldeg; van Spaandonk et al. 2010a) should present a view 
of the disk that does not change much during the orbit. It is also intriguing that GW Lib 
did not show the series of rebrightening events (echo outbursts) following its steep optical 
dechne that are often observed in WZ Sge systems (Patterson et al. 2002). 

A year later, when GW Lib was again observable, the GALEX data begin to show a 
long (hrs) timescale variation that becomes even more pronounced in 2009 and 2010. In 2008 
and 2009, the amphtude in FUV is larger than the NUV, implying the innermost hot areas 
of the disk are involved in the cause of this variation. On 2008 March 29 (Figure 3b), a long 
term modulation that is not well-defined (as it is about the length of the 3 hr light curve) 
appears in the optical data, and there is an even more puzzling feature: a distinct pulse-like 
periodicity. The DFT (Figure 4) indicates a period of 19.08±0.03 min. Copperwheat et 
al. (2009) also report a consistent signal from March 30-April 29 at 19.2±0.2 min with an 
amplitude in g of 2.2%, which declines in amplitude throughout their April data. In addition. 
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Schwieterman et al. (2010) show a similar (19.7-19.9 min) but not identical period in 3 out 
of 8 nights from 2008 May 30 to June 13. As the frequency is not stable, we regard this as a 
quasi-periodic feature that is similar in timescale to those that have been reported in several 
novalike systems (Garnavich & Szkody 1992; Taylor et al. 1999; Andronov 1999). 

After the disappearance of the 80 min and 19 min phenomena, the optical light curves 
in 2008 June and throughout 2009 exhibit a longer modulation near 4-5 hrs (Figures 3c, 3d) 
with amplitude of 0.2-0.3 mag peak-to-peak which is neither consistent nor coherent. This 
timescale is about double the 2.1 hr period first discovered by Woudt & Warner (2002) and 
then studied further in the prc-outburst data of Coppcrwheat et al. (2009). Figure 3d shows 
that while this modulation is very prominent on 2009 May 17, it completely disappears on 
May 18 and 19 (as evidenced on data obtained on two different telescopes). Both GALEX 
and optical data exist on 2009 May 16 (although not simultaneously), showing that a similar 
modulation is apparent in both (Figures 2a and 3c and Table 3). However, while the periods 
are similar on the same date (the combined optical data from May 14-17 fit with a period 
of 191.5±0.2 min and the GALEX data give periods of 218±llmin (FUV) and 210±9 min 
(NUV)), the phases do not seem to match up. The peak of the GALEX data occur near 
730-740 UT, which would imply the previous peak occurred at 530-538 UT, while the optical 
peak is closer to 450 UT. While the origin of this modulation is not understood, the fact 
that it changes phase on timescales of days and sometimes disappears entirely, points to an 
origin in the accretion disk. If it is related to an eccentricity of the disk, the changes mean 
that the disk is still not in a steady configuration two years after the superout burst. 



4. Pulsation Search 

A period search was conducted on all the optical and ultraviolet data, primarily to search 
for the return of the pre-outburst pulsations. For this, the hght curves were converted to 
a fractional amplitude scale by dividing by the mean count rate and then subtracting one. 
Then a Discrete Fourier Transform (DFT) was computed for each light curve up to the 
Nyquist frequency (half the sampling frequency) . The best fit periods and amplitudes were 
found by using linear and non-linear least squares fitting. To find the significance of any 
period, the 3cr white noise limits were determined by using a shuffling procedure (Kepler, 
1993; Mukadam et al. 2010). This involves shuffling the intensities of a light curve (with 
all the periodicities subtracted out), while keeping the times constant in order to produce 
a pure white noise lightcurve. The DFT of this lightcurvc is then computed to find the 
average amplitude between and the Nyquist frequency to yield a la measure. Repeating 
this process 10 times produces an empirical, reliable 3(7 limit. 
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Table 3 lists all the 3a noise levels and any significant results obtained for the ultraviolet 
and optical data. There is no consistent period that is evident as a non-radial pulsation up 
to 3 years past outburst in cither the ultraviolet or the optical data. The amplitudes of the 
650s period evident prior to the outburst were 100 mma in the UV and 10-15 mma in the 
optical. Figures 5 and 6 show the DFTs of the latest data. 

Copperwheat et al. (2009) also obtained photometry in 2007 June and July and 2008 
Mar- June. In their last dataset (2008 June 21), they note a period of 290s with an amplitude 
of 1.25%. There is no consistent evidence for a pulsation near this period in our longest 
datasets (Table 3). Thus, this period is likely associated with flickering noise. 

Van Spaandonk et al. (2010b) recently used VLT spectroscopy to determine a spin 
period of 97±12s from an analysis of a Mg II4481 absorption feature. In contrast to the 
obvious presence of the spin period in the accreting pulsator V455 And (Araujo-Betancor et 
al. 2005), there is no indication of this spin period in any of our data on GW Lib. 

5. White Dwarf Cooling 

The UV data obtained following the outburst of the dwarf nova WZ Sge showed that 
its white dwarf takes more than 3 yrs to cool to its pre-outburst temperature (Godon et al. 
2006). As with all dwarf novae near the peak of outburst, the UV flux is dominated by the 
accretion disk and it is only after the flux has dropped down from the plateau region of the 
outburst that the white dwarf becomes visible. Nogami et al. (2009), Hiroi et al. (2009) and 
van Spaandonk et al. (2010a) discuss the spectra of GW Lib from pre-outburst to the slower 
decline phase up to 3 months post outburst. Our spectrum (Figure 7) obtained one year 
after outburst, is similar to quiescent spectra in showing the absorption from the white dwarf 
flanking the emission from the accretion disk, but there is a steeper blue continuum than at 
quiescence. As the emission lines are relatively narrow and single-peaked, the inclination of 
GW Lib is generally taken to be low (~11°) and the mass to be high (near 1 M©; Szkody et 
al. 2000, Thorstensen et al. 2002, Steeghs et al. 2007, van Spaandonk et al. 2010a,b). 

We modelled the optical spectrum of GWLib using a two-compoment approach. The 
emission from the disk is described with an isothermal and isobaric slab of hydrogen (Gansicke 
et al. 1997). For the white dwarf, we used a grid of model spectra spanning a wide range in 
effective temperatures and surface gravities, computed using TLUSTY195 and SYNSPEC45 
(Hubeny & Lanz 1995) and adopting metal abundances of 0.1 the solar values, as determined 
from our previous HST/STIS observations (Szkody et al. 2002), and using the distance of 
104±2o as determined by Thorstensen (2003). Free parameters in this model are the temper- 
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ature and column density of the hydrogen slab, which determine the equivalent widths of the 
emission lines as well as the Balmer decrement (see Gansicke et al. 1999 and Rodriguez-Gil 
et al. 2005 for more details), and the white dwarf effective temperature and surface gravity. 
The surface gravity, along with the flux scaling factor between models and observation and 
a white dwarf mass-radius relation (Wood, 1995) give the white dwarf mass. 

The disk parameters are adjusted to reproduce the observed H/3 and H7 fluxes, and 
the white dwarf parameters are varied to fit the broad Balmer absorption lines as well as 
the slope of the continuum. As the ionization fraction of hydrogen depends both on the 
temperature and the pressure within the atmosphere, the white dwarf effective temperature 
and its surface gravity are strongly correlated parameters, and acceptable fits are found for 
a 15000-25000 K white dwarf, corresponding to masses of 0.57-0.97 M©. The best fit was 
achieved for a O.97M0, 25,000K white dwarf (Figure 8). Even at this temperature, the 
bluest wavelengths are not fit well, implying some other component may be contributing. 

Although the GALEX and optical data were not simultaneous, we tried to constrain any 
hot component by combining the closest datasets. Figure 9 shows the CTIO data from April 
17 plotted with the GALEX data from May 13 (fiuxes are plotted in frequency to show the 
data better), with the observed range of variabihty of the GALEX daXdi,. The UV data does 
not follow the upturn at the bluest optical wavelengths. A hot 25,000K model lies far above 
the GALEX Y>ovais. This could mean that either there is a large reddening (not very feasible 
given the small distance of GW Lib), or there was a large change in the GALEX or optical 
fiuxes between the two measurements (also not very feasible), or the model parameters are 
not correct. Using the best determination of mass from van Spaandonk et al. (2010) i.e. 
mass of 0.84 M0 and the distance of 104 pc with a small amount of reddening (E(B-V)=0.03) 
yields the fit shown in Figure 9 for a temperature of 19,000K. In this model, the optical sum 
of the white dwarf and disk falls below the observed values. However, this could be possible 
if the disk contributes a greater portion of the optical continuum fiux. The fitting of UV 
spectra have shown that the white dwarfs in accreting pulsators contribute 75-89% of the 
UV light and 42-75% of the optical (Szkody et al. 2010) during quiescence. In the model 
shown in Figure 9, the white dwarf contributes 80% of the fiux near 45OOA. Since GW Lib 
was not yet at quiescence, the disk contribution could be larger than normal. 

Given all the uncertainties (10% in optical fiux calibration, 20-30% in distance, unknown 
changes in the month between the GALEX and CTIO observation times, unknown reddening 
and unknown disk contribution and fiux distribution), it is not worth pursuing a more 
detailed fit at this time. It is clear that there is a hotter component (likely the white dwarf 
as the dominant contributor to the fiux) than quiescence a year after its outburst. In WZ 
Sge, the white dwarf was about 4000K hotter than quiescence one year after its outburst 
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(Godon et al. 2006). Given the similar physical characteristics of the two systems, including 
both having massive white dwarfs, it is interesting to determine if GW Lib remains hotter 
at the same time post outburst than WZ Sge. Since the cooling is related to the amount 
and depth of the accreted mass, the larger outburst amplitude of GW Lib (9 mags) vs WZ 
Sge (7.5 mags) may indicate that more material was accreted during the outburst of GW 
Lib. In addition, the much lower inclination of GW Lib provides a better view of the disk 
radiation at outburst. The determination of the actual cooling of the white dwarf in GW 
Lib compared to WZ Sge will come from ultraviolet spectra. 

A soft X-ray component was seen with Swift (Byckling et al. 2009) within a few days 
of outburst that was consistent with an optically thick boundary layer. The hard X-ray flux 
was abnormally large at outburst (1000 times its pre-out burst value) and was still an order 
of magnitude larger than pre-outburst through 2009 March. Whatever caused the increased 
X-ray luminosity could result in a slower cooling of the white dwarf in GW Lib compared to 
WZ Sge. 

6. Conclusions 

Our UV and optical data on GW Lib during the three years following its outburst has 
shown the following unique traits: 

• The accretion disk shows large optical variability following the steep dechne part of its 
post-outburst light curve. 

• A superhump modulation evident only in optical, not ultraviolet, follows this transition. 

• One year post-outburst, a 19 min quasi-period is visible in the optical for 2-4 months. 

• The ultraviolet begins to show variability on ~4 hour timescales a year after the out- 
burst and the amplitude of this variability increases in the second and third years after 
outburst. 

• The optical begins to show a 4 hr variation 2 years after outburst which is similar to 
the ultraviolet but not in phase and this variation can disappear from one night to the 
next. 

• At one year past outburst, the white dwarf likely remained hotter than at quiescence, 
although a range of temperatures and disk components is possible. 
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• There is no evidence of the return of the pre-outburst pulsations for the 3 years following 
outburst. 

The very large outburst amplitude, the long-lasting X-ray emission, the hot white dwarf, 
the lack of post-outburst rebrightenings, the presence of quasi-periodicities and the long pe- 
riod modulation all point to a large accretion episode, resulting in a hot white dwarf and a 
massive accretion disk. While the structure of the disk that results in all the periodicities 
observed is not known, following the detailed transitions in systems like this provide some 
clues to the types of phenomena that the disk undergoes as it transitions from its dominant 
optically thick emission state to its minor quiescent contribution to the system light. Ob- 
taining the pulsation periods and amplitudes when the white dwarf cools enough to resume 
its modes will reveal the effect of the accretion event on the interior layers of the white dwarf. 

It is interesting to compare the behavior of GW Lib after its outburst to that of two 
other similar systems that also had recent outbursts and contain pulsating white dwarfs. 
SDSS0745-I-45 had an outburst in October 2006, showed no visible pulsations one year post 
outburst (Szkody et al. 2010) but pulsations at the same pre-outburst periods had returned 
by February 2010 (Mukadam et al. 2010b). V455 And had an outburst in September 
2007 and its pulsation return is currently being monitored. A comparison of the outburst 
amplitudes and return timescales and characteristics of the pre-and post outburst pulsations 
of these three systems should provide some clues to the depth of heating and interaction of 
the accretion event with the driving zone of the pulsations. 

This research was funded by NASA GALEX grants NNX08AU43G, NNX09AF87G and 
NSF grant AST0607840. 
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Table 1. Summary of GALEX Data 



UT Date 


mid UT 


Exp {&] 




FUV 




NUV 


2007 May 26 


18:19:55 


1088 


14. 


.439±0.004 


14. 


.874±0.004 




18:58:36 


1120 


14. 


,446±0.004 


14. 


,840±0.004 




21:37:16 


1136 


14. 


,420±0.004 


14. 


,818±0.004 


2007 Jun 08 


13:50:20 


106 


14. 


,70±0.01 


15. 


,07±0.01 


2007 Jun 18 


03:40:31 


221 


14. 


,79±0.01 


15. 


,16±0.01 




05:19:05 


616 


14. 


,799±0.007 


15. 


.165±0.006 




06:57:41 


1002 


14. 


.805±0.006 


15. 


.179±0.005 




08:36:16 


1388 


14. 


.806±0.005 


15. 


.192±0.004 




10:14:53 


1690 


14. 


.780±0.004 


15. 


.180±0.004 


2008 May 13 


17:52:59 


1321 


15. 


.948±0.009 


16. 


.179±0.007 




19:31:39 


1343 


16. 


,046±0.009 


16. 


,222±0.007 




21:12:16 


1201 


15. 


,864±0.008 


16. 


,110±0.006 




22:51:21 


1197 


15. 


,968±0.009 


16. 


,136±0.007 


2008 Jun 03 


11:50:28 


1684 


15. 


.986±0.007 


16. 


.175±0.006 




13:29:03 


1667 


15. 


.917±0.007 


16. 


.142±0.006 




15:07:39 


1678 


16. 


.074±0.008 


16. 


.261±0.006 


2008 Jun 13 


10:05:33 


1683 


16. 


.071±0.008 


16. 


.276±0.006 




11:44:09 


761 


16. 


,02±0.01 


16. 


.21±0.02 


2009 May 16 


10:41:23 


1498 


16. 


,23±0.008 


16. 


,31±0.005 




12:19:59 


1498 


15. 


,88±0.007 


16. 


,03±0.004 




13:58:36 


1503 


16. 


.27±0.008 


16. 


,34±0.005 


2010 Apr 25 


03:35:27 


1488 






16. 


,35±0.005 




05:14:03 


1489 






16. 


.13±0.004 




06:52:40 


1460 






16. 


.50±0.005 




08:35:13 


1471 






16. 


.04±0.004 
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Table 2. Summary of Optical data 



UT 


Site 


JD-2450000 


UT range 




Exp (s) 


Filters 


Comments 


2007 Apr Id 


ATA /TOTT 


A O A/C 

4206 


AO CA AO 1 1 AO 

08:50:03-11:02: 


1 C 

:15 


C A 

60 


T T 
U 






2007 Apr 17 


ATA yrOTT 


A OAT 

4207 


AC CA C1 1A OO 

06:59:51-10:23: 


O/? 

:26 


A 1 01 A C 

4.1-219.5 


T T 
U 






zL)U7 Apr 18 


ATA TOT T 

JNMbU 


A O AO 

4208 


AT A1 AO 1A CO 

07:01:08-10:58: 


:14 


o /I /I r* T 

3.4-45.7 


U BV Hi 


DFT in 


U 


OAAT A 1 A 

zUUy Apr 19 


ATA TOT T 


A onn 
4209 


AT. AO . OA AA. C A 

07:02:o9-09:50: 


. C A 

:o4 


O T AO a 

O.7-9O.0 


T T DT A D T 

U nV HI 


DFT in 


u 


zUU7 Apr zU 


ATA TOT T 
IN Mb U 


/I 01 A 

4210 


CiC . AC A^ 1A./in, 

0d:4o:42-10:49; 


.OO 

:cS2 


O T /I T n 

o. 7-47. 9 


TT DT/" D T 

U JdV HI 


DFT in 


u 


on AT A ^™ oo 

zU(J7 Apr zz 


ATA TOT T 


O 1 o 

4212 


AA AC Ar\ 1A OO 

09:05:40-10:23 


01 

:21 


O C CA T 

3.5-50.7 


U BV Hi 


DFT in 


u 


OA AT A O/^ 

z(J(J7 Apr ZD 


ATA TOT T 

i\l Mb U 


4215 


AT O /I or* AT OO 

07:24:36-07:32 


A A 

:44 


O A OO A 

2.0-33.4 


U BV HI 






OA AT A OT 

zU07 Apr 27 


ATA TOT T 

IN Mb U 


O 1 T 

4217 


r\C; CO CO AT A1 

0o:52:53-07:01: 


or' 

:2o 


O C A A C\ 

2.5-44.0 


U BV HI 






zUU7 May U4 


ATA TOT T 


A OO A 

4224 


ACAA.AO CiCX O, 

Od:09:Oc5-Od:18; 


.AT 

:07 


o c CO n 


U BV HI 






zUU7 May U7 


ATATOTT 


A OOT 

4227 


A C O T. /I 1 AA. A A 

0o:o7:41-09:40: 


. A O 

A6 


1 A 01 A O 

1. 4-219. o 


U BV HI 


DFT in 


u 


OA AT "I\ /T^ 1 A 

zUU7 May ID 


ATA TOT T 


/I oo A 

4230 


AC CA CO A A C /I 

05:50:58-09:54: 


A C 

:4d 


O A 1 O O T 

3.0-138.7 


U BV HI 


DFT in 


u 


O A A T 1\ 1 1 

zUU7 May 11 


ATA TOT T 

In Mb U 


4231 


AO /I O OC A A 1A 

08:42:35-09:19 


C A 

:50 


O 1 OCA 

3.1-85.0 


U BV HI 


DFT in 


u 


O A A T A T„ 1 O 

2007 May lo 


ATA TOT T 

JNMbU 


4233 


AO 1 T AO A A OO 

08:17:08-09:33: 


1 C 

:15 


/( O 1 £?0 T 

4.3-162.7 


T T DT Z T> T 

U BV Hi 


DFT in 


u 


OA AT A /T^ 1 /I 

zUU7 May 14 


ATA TOT T 

INlVlbU 


/I OO /I 

42o4 


AT CO CO A A OC 

07:52:53-09:25 


1 O 

:18 


O O 1 AT A 


T T ryx z T) T 

U BV HI 


DFT in 


u 


OA AT A/T^ OA 

2007 May 29 


TX TTA/'AT 

WlYiN 


A^ Af\ 

4249 


AC AO AA AT CT 

05:43:44-07:57: 


:11 


t c 

15 


T> O A 






2007 Jun 01 


ArU 


/( O CO 

4252 


AO A A CC A A AO 

08:09:56-09:02: 


c c 

:56 


1 A 

10 


T)f~^ ACi 






2007 Jun 05 


ATA TOT T 

JNMbU 


4256 


AC A1 AC AC 1 O 

05:01:46-05:13: 


OA 

:30 


A A AO O 

9.4-92.8 


T T DT 7" T> T 

U BV HI 






OA AT T r\l2 

2007 Jun Oo 


ATA TOT T 

IN Mb U 


,1 O C T 

4257 


r\ A CC O /I A/^ Of 

04:55:34-0o:3d: 


:oO 


01 OOO /I 

21.D-228.4 


T T DT D r 

U BV HI 






OAAT 1 , , , , AO 

2007 Jun Oo 


ATA TOT T 

iNMbU 


/I o c n 
4259 


AO.OO./11 AC. OA 

0o:o2:41-05:20: 


. A C 

:4o 


/I A A1 /I 

4.4-91.4 


T T DT A D T 
U BV HI 






z007 Jun 14 


ATATOTT 
INMbU 


A OC C 

42o5 


AO. /IT. AT C\CA'\. 

0c5:47:07-0d:41; 


. A n 
:49 


O O 1 c c c 

8.2-155.5 


TT DT r T~>T 

U BV HI 






OAAT T 1 C 

2007 Jun 15 


ATA TOT T 

IN Mb U 


42dd 


AO CA OO AT 1 O 

03:59:38-07:13 


c o 

:52 


O 1 1 O /I O 

0.1-124.8 


T T DT /" O T 

U BV HI 






OAAT T 1 O 

2007 Jun 18 


ATO A TT T 

iNbMU 


4259 


AC OA AC /1A 

05:30:39-05:40 


o o 

:28 


T C O 1 

7.5-63.1 


T T DT /" D 7" 

U BV HI 






O A A T T 1 A 

2007 Jun 19 


ATA TOT T 

IN Mb U 


A OTA 

4270 


A /I CA O*^ AC AA 

04:59:26-05:09 


C A 

:59 


5.3-103.0 


U BV HI 






200» Mar 29 


Ar^U 


4oo4 


nT.cc/io 1 1 .OC 

07:56:43-11:25: 


:o3 


10 


r$(j40 






2008 Apr 17 


CTIO 


4573 


07:08:19-08:09: 


:19 


2x1800 




spectra 




2008 Jun 13 


NMSU 


4630 


03:51:41-07:03: 


:03 


6.7-136.6 


UBVRI 


DFT in 


u 


2008 Jun 20 


NMSU 


4637 


04:01:36-06:16: 


:10 


3.0-112.3 


BG40,UBVRI 


DFT in 


BG40 


2009 May 14 


NMSU 


4965 


23:21:59-01:16: 


:45 


60 


BG40 






2009 May 15 


NOFS 


4966 


05:15:10-09:14: 


:30 


20 








2009 May 16 


NMSU 


4967 


01:11:23-03:24: 


:12 


60 


BG40 






2009 May 16 


SRO 


4967 


04:34:30-09:50: 


:26 


45 


clear 






2009 May 16 


NOFS 


4967 


06:04:55-09:45: 


:49 


45 


B 






2009 May 17 


NOFS 


4968 


04:51:22-09:19: 


:45 


45 


B 






2009 May 17 


SOAR 


4968 


05:04:24-09:58: 


:28 


8.9 


V 






2009 May 18 


SAAO 


4969 


21:47:40-01:39: 


:30 


40 


B 






2009 May 19 


NMSU 


4970 


05:19:52-08:28: 


:29 


60 


BG40 






2010 Mar 04 


KPNO 


5259 


10:07:09-12:36: 


:21 


15 


BG39 
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Table 2 — Continued 



UT 


Site 


JD-2450000 


UT range 


Exp (s) 


Filters Comments 


2010 Mar 05 


KPNO 


5260 


10:03:19-11:36:37 


15 


BG39 
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Tabic 3. STimmary of Observed Periods and Limits 



TTT 
U ± 


uata 


TM 


3(7 (inma) 


Periods (min) 


Amp (mma) 


070416 


NMSU 


28 


14.0 






070417 


NMSU 


145 


9.2 






070418 


NMSU 


199 


6.9 






070419 


NMSU 


70 


11.2 






070420 


NMSU 


192 


5.6 






070422 


NMSU 


52 


3.6 






070510 


NSMU 


96 


49.1 


80.1±1.3 


100±9 


070514 


NMSU 


41 


73.9 


83.3±1.3 


137±3 


070526 


GALEX NUV 


120 


17.1 


32.3±0.4 


21±3 


070526 


GALEX FUV 


120 


26.9 






070529 


WIYN 


311 


12.8 


82.9±0.9 


52.8±1.5 


070601 


APO 


278 


14.3 






070618 


GALEX NUV 


294 


36.2 


154±2 


78±7 


070618 


GALEX FUV 


294 


44.4 






080329 


APO 


1255 


5.7 


204±4, 85.2±1.8,19.08±0.03 


39.5±0.8,14.7±0.7,24.9±0.8 


080513 


GALEX NUV 


88 


33.6 


327±9 


67±6 


080513 


GALEX FUV 


88 


52.5 


291±9 


96±10 


080603 


GALEX NUV 


87 


39.3 


275±22 


66±7 


080603 


GALEX FUV 


87 


49.5 


219±8 


106±22 


080613 


GALEX NUV 


42 


38.7 


215±27 


55±18 


080613 


GALEX FUV 


42 


48.6 






080620 


NMSU 


89 


27.1 


124±6 


53±6 


090514-15 


NMSU 


78 


34.3 






090515 


NOFS 


719 


10.9 


248±7 


60±2 


090516 


NMSU 


90 


47.7 


308±108 


144±50 


090516 


SRO 


187 


49.2 


293±10 


128±5 


090516 


NOFS 


169 


41.2 


206±11 


109±7 


090516 


GALEX NUV 


78 


81.7 


210±9 


150±8 


090516 


GALEX FUV 


78 


105.7 


218±11 


191±9 


090517 


NOFS 


200 


27.6 


283±30 


78±11 


090514-17 


multi 


1438 




192±3 


54±2 


090518-19 


SAAO 


239 


8.0 






090519 


NMSU 


61 


20.5 






100304 


KPNO 


132 


13.9 


79.0±3.5, 40.2±0.9 


21±3, 23±3 


100305 


KPNO 


90 


16.2 






100425 


GALEX NUV 


30.5 


103 


229±4 


205±8 
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Fig. 2a.— GALEX NUV (left) and FUV (right) light curves. Each point in 2007 data is 
a 29s integration, in later data is a 59s integration with error bars for each orbit listed in 
Table 1 . The magnitudes are on the AB system as explained in the text. 
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Fig. 2b.— GALEX NUV light curve in 2010. Details as in 2a. 
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Fig. 3a. — Individual optical light curves obtained from 2007-2010. Error bars are generally 
a few hundredths of a mag. Filters are as listed in Table 2, with filter same as for DFT if 
more than one on a night. Dates increase across each row and then down the column. 



-24- 



o 







" a 




- 2 




- 1^ 




- o 




o 




- 1:^ 








- Ln 




- 2 




z 









U I I I I I 



O 

o 

CD 



O 

in 



o 
o 

lO 



o 
o 



o 
- in 



CO 




o 



o 

o V- 



o 
in 



-co 
o 
o 



o 
in 
in 



o 
o 
in 



O 3 

in c 



o 
o 



o 
in 



1^ oD cn 
in in in 



CD CD 



in 

CD 



CD 
CD 



1^ 
CD 



sprniu5D|A| 



spn:i]u6D|A| 



sprn]u6D|/\| 



1^ 
o 
o 

(N 

(/I 



V 



o 
in 



o 
o 
in 



in 



o 
o 



o 
in 



o 
o 

-CM 



CD 



1^ 



00 



o 

CM 
CO 



o 
o 

CO 



o 

CO 

in ■ 



o 

CD 

in 



o 
in 




aprniu5DiA| 



8pn:)iu6D|A| 



8prniu6o|/\| 



Fig. 3b. — Continued. 
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Fig. 3c. — Continued. 
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Fig. 3d. — Continued. 
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Fig. 4. — Fractional intensity (top), DFT (middle) and enlargement of low frequencies (bot- 
tom) of APO optical data on 2008 March 29 showing the 19.08 min (1145s) period. Dashed 
line is the 3a noise limit. 
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Fig. 5. — Fractional intensity (top), DFT (middle) and DFT with 3.822 hr period removed 
(bottom) of 2010 April 25 GALEX NUV data. Dashed lines show the Scr limits for the noise. 



-29- 




I I I I I I I I I I I I I I I I I I I I K'^Ti |yyvi I I I I I I I I I I 

0.001 0.002 0.003 0.004 0.005 0.006 0.007 

Frequency (Hz) 



Fig. 6. — Fractional intensity (top), DFT (middle) and window function (bottom) of 2010 
March 4 optical data along with the 3a noise limit (dashed line). The pre-outburst pulsations 
were all at frequencies larger than 0.001 Hz. 
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Fig. 7. — CTIO spectrum of GW Lib. Note the broad absorption from the white dwarf 
flanking the emission from the accretion disk. 
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Fig. 8. — CTIO spectrum (grey) fit with a 25,000K, log g=8.55 white dwarf at 104 pc (red) 
and an optically thin emission line component (green). 
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Fig. 9. — GALEX data range from 2008 May 13 (solid points) combined with the CTIO 
spectrum from April 17 (grey) fit with a 19,000K, log g=8.36 white dwarf at 104 pc with a 
reddening of E(B-V)=0.03 (red line) and an emission line component (green). 



